Noise Properties of Coherent Perfect Absorbers and Critically-coupled Resonators 
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The performance of a coherent perfect absorber (time-reversed laser) is hmited by quantum and 
thermal noise. At zero temperature, the quantum shot noise dominates the signal for frequencies 
close to the resonance frequency, and both vanish exactly at the resonance frequency. We compute 
the sensitivity of the absorbing cavity as a background- free detector, limited by finite signal or 
detector bandwidth. 



In recent work the authors and collaborators have pro- 
posed [ij and demonstrated [2 the phenomenon of coher- 
ent perfect absorption, or "time-reversed lasing". Ap- 
plying the time-reversal operation to the classical elec- 
tromagnetic equations yields the following statement: If 
a cavity containing a gain medium characterized by an 
amplifying refractive index n(r) = ni(r) — m2(r) reaches 
the lasing threshold for a certain lasing mode at fre- 
quency cjo, then a cavity containing a dissipative refrac- 
tive index n*(r) will perfectly absorb an input mode at 
ujQ corresponding to the time-reverse of the lasing mode. 
We refer to such a lossy cavity as a coherent perfect 
absorber (CPA). Assuming the input signal is perfectly 
monochromatic, the above "CPA theorem" is rigorously 
true within classical electromagnetic theory, where the 
effects of quantum and thermal noise are neglected. The 
CPA is a generalization to arbitrary geometries and arbi- 
trary numbers of input channels of the well-known con- 
cept of a critically-coupled resonator (CCR) [3 , an opti- 
cal device which can be used for switching, modulation, 
enhanced photodetection and sensing [H |5] , and which 
may be regarded as the single channel limit of a CPA. Be- 
cause a CPA is associated with a vanishing output signal 
in the classical zero-temperature limit, a CPA (or CCR) 
can function as a background-free detector or interferom- 
eter, similar to a Mach-Zender interferometer (MZ). The 
fundamental limits to its effectiveness in this role are de- 
termined by quantum and thermal noise, which are the 
subject of the present paper. 

Quantum fluctuations, such as spontaneous emission, 
break the symmetry between emission and absorption. 
Hence, noise processes in a CPA (or CCR) differ from 
those in a laser or amplifier. Within semiclassical theory, 
the laser has zero linewidth; by including quantum fluctu- 
ations, one obtains the Schawlow-Townes (ST) linewidth 
^ , which decreases inversely with the output power well 
above threshold. As is well-known, the ST linewidth 
arises from the dephasing of the above-threshold laser 
field due to quantum noise (usually characterized as "one 
noise photon per relaxation time per mode"). Because a 
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CPA does not contain an inverted medium like a laser, 
and spontaneous emission vanishes at low temperature, 
there is no direct analog of the ST linewidth in a CPA. 
At zero temperature, the frequency characteristics of the 
output are determined entirely by the input, and the 
CPA absorption resonance has approximately twice the 
passive cavity linewidth, as expected for critical coupling. 
However, this linewidth is further modified, as we shall 
show below, by an analog of the Petermann factor [7l[8]; 
this effect has not, to our knowledge, been recognized in 
absorbing systems. 

Due to the lack of spontaneous emission at zero tem- 
perature, were it possible to probe the resonance with 
a truly monochromatic input at the CPA resonance, the 
perfect absorption predicted by the classical theory would 
occur — even with quantum fluctutations taken into ac- 
count. In this situation, the quantum shot noise van- 
ishes simultaneously with the zero of the average output. 
In practice, however, the unavoidable linewidth of the 
input field, Ain, combines with quantum shot noise to 
generate a finite noise floor even at the resonance fre- 
quency. The noise dom inates the signal within an in- 
terval 8uj^ oc Ain /Pin around the resonance frequency, 
where Pin is the power of the input signal. This be- 
havior is similar to a MZ, with the absorption into an 
external reservoir playing the role of an unobserved MZ 
output port. However, in contrast to the MZ, the res- 
onant frequency response of the cavity can give rise to 
parametrically better signal-to-noise ratio in the vicinity 
of the background-free point [9 . For finite temperature, 
the output also includes thermal noise; at the absorption 
resonance, the thermal emissivity takes the black-body 
value for a one-port system (CCR), but can be signifi- 
cantly reduced in a multi-port CPA. 

We can analyze the noise properties of the CPA using 
the standard input-output framework of quantum optics 
[10, 11 . Input and output photon operators are con- 
nected by the classical electromagnetic scattering matrix 
Sij{uj)^ where j = 1, 2, . . . , A/" denote the scattering chan- 
nels. The ^-matrix for a CPA is sub-unitary due to 
the presence of an absorbing reservoir, and describes a 
"gray-body" , with some fraction of the incident photons 
absorbed, and the remainder scattered out at the input 
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frequency. The scattering channels denote different spa- 
tial states in the asymptotic region which suffice to rep- 
resent an arbitrary incoming or outgoing field at cj, e.g. 
incoming and outgoing angular momentum channels in a 
two-dimensional scattering geometry. The input photon 
operators are denoted by aj , and the output photon oper- 
ators by bj . The special case of a single scattering channel 
{N = 1) corresponds to a standard CCR. We assume that 
the CPA is coupled to an ideal external reservoir, so that 
the absorption of light produces negligible heating and 
hence negligible change in the output noise. The photon 
operators are related by an input/output relation [ 101412] 

bi{uj) = ^Sij{uj)aj{uj) ^^Uij,{uj)cj,{uj), (1) 



where the c^^'s are ladder operators for reservoir quanta. 
Here we have assumed a coupling to the reservoir which 
adds the minimum amount of quantum noise [13 . The 
requirement that a, 6, and c obey canonical commutation 
relations, e.g. [ai{uj)^ a^j{uj')] = Sij S{uj — uj')^ yields the 
fiuctuation-dissipation relation [TT] 



SS^ ^UU^ = 1, 



(2) 



where 1 is the N x N identity matrix. Equation ([2| gen- 
eralizes the unitarity relation of the lossless system, and 
implies that the eigenvalues of the S'-matrix generically 
have magnitude smaller than unity. 

We are interested in the shot noise in the output field, 
for a coherent input at some frequency uj. To work with 
equal- frequency correlators, it is convenient to rescale 
the continuum operators {ai^bi^Cjy} to discrete opera- 
tors {di^bi^Cj^}^ which are normalized so that the equal- 
frequency commutator is unity, e.g. [di{uj), aj{uj)] = Sij. 
Here and in the following, we omit the uj dependence from 
the notation. Next, the input photon operator can be re- 
written using a displacement transformation [131 US 

= + a-, (3) 

where ai is a coherent state amplitude, and a- is an op- 
erator accounting for fiuctuations around the coherent 
state, which likewise obeys the canonical commutation 
relation [a^, ah = 



Hence, 



(4) 



The operator A/^ = blbi gives the output photon fiux 
per unit frequency (at frequency uj) in channel i. (Cor- 
respondingly, Xi = huJoJVi describes the spectral density.) 
Using Q, we can calculate the expectation value and 
correlation function for A/^, using standard Gaussian n- 
point operator correlators [15]. We take {d[) = (cj^) = 0, 
and (a^ctj) = (zero net fiuctuation around the specified 



/(T) = [ex.p{huj/kBT) — 1] and T is the temperature 
of the reservoir. The result for (Afi) is 

{M) = \{Sa),\'+[l-{SS^)u] f{T). (5) 
Thus the total output is 

i 

and the noise is 
(5AA2)^^[(MAG-)-(M)(AG-)] 

= \Saf 
+ f{T){2\Saf -2\S^Saf 

+ [fiT)fTv\{l-SS^f 



Tr 



(1 - SS^) } 



(7) 



For T ^0, f{T) 
nian result 



0, ^ and (|7| reduce to the Poisso- 



{N) = {5N^) = \Sa\' 



(8) 



The CPA condition is achieved when the index of refrac- 
tion of the cavity, n(r), is chosen such that there exists 
an eigenvector of the ^-matrix with eigenvalue zero at a 
specific input frequency ujq. If a is chosen to be this eigen- 
vector, then we see from Eq. ([8| that the mean outgoing 
photon fiux vanishes (as it should be from the classical 
CPA theorem), and so does its variance (the shot noise). 
This confirms the statement that the CPA effect is unaf- 
fected by quantum fluctuations at zero temperature for 
a purely monochromatic input. 

At finite temperature, the shot noise dominated regime 
corresponds to T <C hujo/ks, For T 7^ 0, (AT) and (W^) 
are unequal, and differ by a term arising from the beating 
between the thermal emission and the scattered fiux: 

(SAf^) - (AT) = /(T) [2\Saf - 2\S^ Saf] + 0{f). (9) 

This difference can be shown to be strictly positive for 
any input field. The term linear in /(T) vanishes when 
the input field, a, corresponds to the zero eigenvector, 
so thermal noise in the CPA is very small in the low 
temperature regime. 

To estimate the sensitivity of a CPA as a detector, 
we study its behavior for frequencies near the perfect 
absorption resonance, ujq. For now, we consider T = 0. 
In the absence of loss, the S'-matrix has a sequence of 
resonances, associated with poles in the lower half of the 
complex UJ plane and symmetrically placed zeros in the 
upper half plane [1 , at discrete frequencies uj = ujq ^ 
ijc/"^' Near cjq, one of the S'-matrix eigenvalues takes 
the following approximate form [17 : 



coherent input amplitude), and (cj^Cjy) 



S^,f{T) where 



s{uj) 



UJ - ujQ- i^cl"^ 



(10) 
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where ip is an irrelevant phase factor and 7c is the cavity 
hnewidth (FWHM). Eq. (10) is a "single resonance ap- 
proximation" which ignores the presence of other nearby 
poles and zeros. It satisfies the requirements that |5| = 1 
when u; G M (the lossless /S-matrix has unimodular 
eigenvalues), and that s goes to zero and infinity at 
u; = cjo i ^7c/2. Adding absorption pushes the zero 
and pole down in the complex frequency plane, as shown 
in Fig. [ija); to lowest order, they move down by equal 
amounts. To achieve the CPA condition, exactly enough 
absorption is added to push the zero down to the real uj 
axis. Then the eigenvalue of the absorbing cavity is 



(11) 



where the parameter 5^ <C 7c represents a small detun- 
ing of the material loss from the perfect absorption res- 
onance. Note that the resonance of the absorbing cav- 
ity has twice the width of the passive cavity resonance 
(7c/2 7c). This is due to the "critical coupling" con- 
dition that the absorption loss rate is equal to the scat- 
tering loss rate. 

In this case, if the input mode, a, corresponds to the 
zero eigenvector of the ^'-matrix, then (at T = 0): 

{M) = (W^) « (;-"°\^/t (12) 

(cJ-CJo)^+7c ^^Ain 

Here we have expressed the input photon fiux per unit 
frequency, in terms of the input power 



Pin(a;) = |a(a;)|^na;Ai3 



(13) 



This power is averaged over some bandwidth Ain, which 
is taken to be the frequency resolution of either the input 
state or the output detector, whichever is smaller. 

According to Eq. (12), as 5^ and uj ^ ujq^ the 



mean and variance of the photon flux behaves as 



(AT) = (W) 



(14) 



The CPA resonance has a quadratic zero at ujq (Fig. [2| , 
and, in the single resonance approximation, its width is 
7c, the critically-coupled cavity linewidth. (We assume 
that 7c is much less than the free spectral range of the 
resonator). 

However, adding a lossy medium to the resonator 
affects the ^'-matrix eigenvalue beyond the single- 
resonance approximation, and the corrections to this ap- 
proximation generically increase the width. This is the 
exact analog of the Petermann factor in lasers, which in- 
creases the ST linewidth [71 [8] . This is most easily seen 
in the one-port case, where we can express the ^-matrix 
eigenvalue with a zero on the real axis as 



n 



(15) 
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FIG. 1. (color online) Scattering properties of a one-port 
coherent perfect absorber (critically-coupled resonator) con- 
sisting of a one-dimensional uniform dielectric slab of length 
L, with a perfect mirror on one side and vacuum on the other. 

(a) Location of poles (squares) and zeros (circles) in the com- 
plex frequency plane. Open symbols show the poles and ze- 
ros for the passive cavity with refractive index n = 3, and 
filled symbols for an absorbing cavity with refractive index 
n = 3 + 0.005i. The dashed and dotted lines are guides to the 
eye for the geometric interpretation of the Petermann factor 
given in Eq. (15). Without absorption, Icj — cj^l = \uj — 
where cj G M and cj^ and cj^ are the frequencies of a neigh- 
boring pair of zeros and poles (green dots). With very nar- 
row band absorption, near a CPA resonance (in this case 
cj ~ 70.1c/L), this remains approximately true. However, 
for broadband absorption, \(jO — (jJ^\ < \oj — (jU^\ (blue dashes). 

(b) Plot of |s(cj)|^, where s{u) is the eigenvalue of the scat- 
tering matrix (i.e. in this one-port case the reflection coeffi- 
cient). The red circles show exact numeric results obtained by 
the t ran sfer matrix method, while the solid blue curve shows 
Eq. (15), with the product taken over the 20 pairs of poles 
and zeros nearest to cj = 70L/c. The dashed curve shows the 



single-resonance approximation, Eq. (11), using the pole and 
zero nearest to = 70L/c. 



where n indexes the other pairs of zeros and poles, which 
have frequencies and respectively. This ansatz 
agrees very well with exact numerical calculations of 
s{uj), as shown in Fig. [ijb). Due to time-reversal sym- 
metry, in the absence of loss, each factor in the prod- 
uct would be unimodular, since the zeros and poles are 
symmetrically distributed around the real axis. The in- 
troduction of loss breaks this symmetry for all the poles 
and zeros, not just the zero at ujq; this causes the ze- 
ros to move towards the real axis and the poles away, 
so that \uj — uj^\ < \uj — uj^\ for any real a;, as shown in 
Fig.jlja). This effect is neglected in the single resonance 
approximation of Eq. Including it leads to 



(AT) = (W> 



{uJ-UJo)^ Pin(^) 



(16) 
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FIG. 2. (color online) Plot mean output power (A/") and shot 
noise power (v^^A/^), as a function of frequency near a CPA 
zero, for T = 0. The cavity is a two-sided uniform dielec- 
tric slab with refractive index 3 + 0.005i and length L (see 
schematic). For this case the CPA eigenmodes are simply in- 
coming coherent plane waves of equal amplitude from each 
direction, with either even or odd parity. Results for an even 
CPA mode at frequency cj ^ 138.23 c/L are shown. We have 
chosen the input intensity = 25. The solid curves are ob- 
tained from Eq. ([s]), using the transfer matrix method to find 
S{(jj). The dashed curves are obtained by averaging these val- 
ues over a bandwidth Ain = 0.1 c/L, corresponding to a finite 
spectrometer resolution. Inset: schematic of the system. 



where K = Yin \ ~ ^n)/(^ ~ ^n)l ^ ^ ^^r the one- 
channel case. Calculating the Petermann factor for the 
time-reversed counterpart of this cavity, which contains 
gain instead of loss and sits at the first lasing threshold, 
would yield exactly the same value of K [TF . This equiv- 
alence follows from the properties of the iS-matrix under 
time-reversal, as does the CPA effect itself. For a more 
complex, multichannel cavity, the explicit calculation of 
the S'-matrix eigenvalues is more complicated, but the 
symmetry that leads to the factor 1/K in the eigenvalue 
still holds. In other words, adding loss to the cavity so 
as to reach the CPA condition increases the critically- 
coupled cavity resonance linewidth from 7c V^7c- 



In calculating this broadening effect, one could assume 
that the loss of the cavity is broadband, affecting all rel- 
evant poles and zeros equally; a similar assumption of 
broadband gain is made in standard calculations of the 
Petermann factor. When this is not the case, both for 
CPAs and lasers the Petermann correction is reduced, 
and needs to be calculated using the frequency-dependent 
refractive index of the cavity [18 . For the subsequent 
analysis, we will assume that the K occurring in Eq. (16) 
is a given parameter. 



Exactly at cjo, both signal and noise vanish; near cjq, 
the shot noise [(JA/"^)]^/^ dominates over the output sig- 
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FIG. 3. (color online) (a) Schematic of a one-channel absorb- 
ing cavity, where the output amplitude 6 is a superposition 
of the input amplitude a and the reservoir operator c. (b) 
The corresponding Mach-Zender interferometer, with a and c 
entering in the two input ports and b being one of the outputs. 



nal (A/"). The crossover frequency scale is 



\a{uJo)\ 



1/2 



vTtc (17) 



The measured values of the signal and noise must, how- 
ever, be obtained by averaging (16) over the bandwidth 
Ain at each frequency, as indicated by the dashed curves 
m Fig. (2] These averaged values do not vanish at = c^o . 
Up to a factor of order unity depending on the averaging 
procedure, their residual values are 



Pin(c4;o) Ain 



(18) 



In particular, we can regard [(^^A/"^)^^]^/^ as the effective 
shot noise level at the absorption resonance. It dominates 
over the bandwidth-averaged signal if Ain ^ ^^x- 

The crossover frequency scale Su^ is related to the sen- 
sitivity of the output to the loss detuning parameter ^7. 
From Eq. (12), the change in the output signal ai uj = ujq 



resulting from J7 7^ is 



(AM) 



(19) 



(To lowest order, this quantity is unaffected by 
bandwidth- averaging.) The sensitivity of the CPA as a 
low-background detector is given by the minimum ^7 for 
which (AA/") is distinguishable from the effective noise 
level. Comparing ([19]) to (p!8|), we obtain the result 



1^71 > ^/SuJ^Ai, 



(20) 



Since a CPA can serve as an absorbing interferome- 
ter/detector, it is useful to compare it to a lossless inter- 
ferometer such as a Mach-Zender interferometer (MZ). 
As shown in Fig. [3j an analogy can be made between a 
CCR (single-channel CPA) and a MZ. The input pho- 
ton operator a and the reservoir operator c, from the 
one-channel version of Eq. ([T]) , map onto two input pho- 
ton operators for the MZ; meanwhile the output pho- 
ton operator b maps onto one of the MZ outputs. The 
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fluctuation-dissipation relation, Eq. ([2|, is equivalent to 
the relation for transmission into the observed port of the 
MZ, 5*5, (which follows from the unitary scattering ma- 
trix for both ports): 



Sb 



(e*^-l), \U'\ = l-\Sh\' 



(21) 



where is the phase difference between the two arms 
of the interferometer. Perfect absorption in the cavity 
corresponds to ^ = 0, so that the a input is directed en- 
tirely into the second, unmonitored output port, leaving 
the b port empty. Quantum fluctuations in the empty 
c input contribute to the shot noise. Note that in the 
CCR, the single input beam interferes with itself and it 
does not function as an optical interferometer, i.e. it does 
not measure the relative phase of two input beams. The 
multi-channel CPA does, however, act as an interferome- 
ter. In the simplest case of two input channels and a cav- 
ity with parity symmetry, the two /S-matrix eigenmodes 
have even and odd symmetry respectively, and only one 
of these are perfectly absorbed at a given absorption reso- 
nance. The interferometry is performed by making small 
phase changes between the two input beams, moving the 
system slightly away from the eigenchannel. For N > 2 
and/or absent parity symmetry, both the magnitude and 
phase of each input amplitude must be tuned in order 
to reach perfect absorption; this is analogous to a multi- 
input MZ-like interferometer whose input amplitudes and 
phases may be tuned to send all the output into a single 
port, playing the role of the absorbing channel. 

The main difference between the CPA and the MZ lies 
in the frequency dependence. For the latter, Sb varies 
sinusoidally in frequency with free spectral range ^ c/L^ 
where L is the dimension of the system. Absorbing cav- 
ities, however, are described by Eq. (11), and for high- 



Q cavities the absorption resonances are much narrower 
than c/L. Since the signal-to-noise ratio depends in- 
versely on 7c, 1/1/ in the two cases, for a given value of 
cj — cjo, the CPA would have a better signal to noise ratio. 
This resonant enhancement of the sensitivity of a CCR 
has been pointed out in a different context in Ref . [9 . A 
similar effect could be achieved in a MZ by adding res- 
onant cavities along each arm to effectively increase the 
optical path length. 

Finally, we make some comments about the T > case. 
In Eq. (6|, we see that {Af) is written as the sum of the 
classical scattered flux and the gray-body thermal emis- 
sion, /(T) N[l — a], a being the mean scattering strength 



per channel [TT]. For a CCR (one-channel CPA), a = 
at the operating frequency of the absorption resonance, 
so the thermal emission has the black-body value. In the 
MZ analogy. Fig. [sj^b), this is equivalent to connecting 
the c input port to a black-body source. Since the MZ is 
tuned so that the a input is completely directed into the 
unmonitored output port, the black-body emission into 
c is directed into the monitored output port b. 

For a multi-channel CPA, the thermal emissivity is 
less than the black-body value, even at the operating 
frequency of the absorption resonance. This is because 
only one of the N scattering strengths vanishes; the other 

— 1 scattering strengths are nonzero, so that a > 0. 
From Eq. ([9|, the same is true of the thermal contribu- 
tion to the output noise (A/"^). Furthermore, for A/" 1 
a "hidden black" scenario is possible, in which the CPA 
perfectly absorbs the input fleld, a, but emits a negligible 
amount of thermal radiation [16 . More speciflcally, for 
a weakly absorbing system it is possible that the mean 
albedo (reflectivity) can be large, a ^ 1, indicating an 
almost white body, and implying that the thermal con- 
tributions to (JV) and (A/"^) vanish, while nonetheless the 
system is perfectly absorbing (down to the quantum noise 
floor) if the correct input fleld is supplied. Elsewhere, two 
of the authors have shown that this effect can be general- 
ized beyond the case of a perfect CPA at resonance, to a 
disordered scattering medium with weak absorption over 
a large range of frequency and absorptivity [16^. 

The presence of thermal noise also affects the sensitiv- 
ity of the CPA as a detector. From Eq. ([7|), the noise 
level at a; = cjo in the thermal noise dominated limit is 



■2 \ thermal 



g{T) « /(T) Tr (l - SS^) . (22) 



Comparing this to the signal (19), we obtain the sensi- 
tivity limit 



(23) 



The crossover between the bandwidth-dominated and 
thermal-dominated noise regimes occurs at 



5(T) 



(24) 
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